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The big question mark: More than four decades after the discovery of pulsars, the 
composition of matter at their cores is still a mystery. This white paper summarizes how 
recent high-precision measurements of millisecond pulsar masses have introduced new exper- 
imental constraints on the properties of super-dense matter, and how continued timing of 
intriguing new objects, coupled with radio telescope surveys to discover more pulsars , migh t 
introduce significantly more stringent constraints. (Figure from iLattimer fc Prakasljl2004l ). 
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Central Questions 



1. What is the nature of matter at supra- nuclear densities? 

2. What is the equation of state for matter in the centers of neutron stars? 

3. What are the upper mass hmit, upper spin hmit, and typical radii of neutron stars? 

4. What is the distribution of masses for millisecond pulsars? 

2 Studying the Properties of Dense Nuclear Matter 

Given their small size {R ~ 10 km) and large mass (M ~ IAMq), neutron stars contain 
some of the densest matter in the Universe. They are unique astrophysical laboratories for 
testing theories of nuclear matter under high pressure. Neutron stars may exhibit condi- 
tions and phenomena not observed elsewhere, such as hyperon-dominated matter, decon- 
fined quark matter, superfluidity and superconductivity with critical temperatures near 10^° 
kelvin, opaqueness to neutrinos, and magnetic fields in excess of 10^^ Gauss. 

The microscopic behavior of nuclear matter, encapsulated in the equation of state (EOS), 
must be inferred from potentially observable, macroscopic neutron star properties, such as 
their masses and radii. This connection can be made by numerically solving the relativistic 
equations of hydrostatic equilibrium (the Tolman-Oppenheimer-Volkov equations) to calcu- 
late mass-radius [M—R) relations for a given EOS. As Fig.[T]shows, there are large variations 
in predicted radii and maximum masses for different candidate EOSs. These reflect basic 
uncertainties about the behavior and composition of matter at and above nuclear density. 
Some EOSs, such as GSl, assume that large percentages of matter are in exotic states (hyper- 
ons, deconfined quark matter), any of which produces a decrease in the nucleonic (protonic 
and neutronic) degeneracy pressure at any given density, since a larger variety of particles 
are present. This means that matter is relatively more compressible; this produces smaller 
maximum stellar masses. Other EOSs, such as MSO, assume a larger fraction of nucleons 
and have higher pressures for any given density. This means that they predict matter to be 
relatively incompressible and result in larger neutron star mass limits. 

3 Constraints on Pulsar Masses and Radii 

There are several ways in which observations of radio pulsars constrain the properties of 
dense matter. (Considerable effort is also being put into observations of X-ray pulsars to 
determine neutron star radii and masses. This is the subject of a separate white paper. X- 
ray instrumentation can also study neutron star cooling, which is important to understand 
nuclear processes and the composition of dense neutron matter well within the stellar core.) 

3.1 Constraints from pulsar spin rates 

Pulsars must be compact so that they do not shed mass due to centrifugal forces. For a pulsar 
with given rotation period and mass, this puts an upper limit on the radius. The strictest 
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Figure 1: Mass-Radius relation for neutron stars. Each black curve represents a family of neu- 
tron stars masses and radii according to a given equation of state. The region bounded by the 
Schwarzschild condition R < 2GM/c^ is excluded by general relativity, and that bounded by 
R < SGM/c^ (labeled "causality") is excluded by requiring the speed of sound inside the star to 
be smaller than the speed of light. The mass-shedding limit for the fastest spinning radio pulsar 
(PSR J1748— 2446ad, with spin frequency 716 Hz) is labeled "rotation"; points in the green region 
below this line are not allowed for that particular pulsar. Stricter constraints may arise from X-ray 
sources like XTE J1739— 285 (dashed curves calculated under different neutron star models) if their 
spin frequencies are confirmed to be higher than that of PSR J748— 2446ad, potentially excluding 
some equations of state (such as GM3) which lie almost entirely below the rotation curve. A recent, 
precise millisecond pulsar (MSP) mass me asurement (for PSR J1903 -I-0327) excludes the "softest" 
EOSs (red horizontal line) (Adapted from iLattimer &: PrakashI 120071 )]. 
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constraint comes from the fastest-spinning p ulsar. The present confirmed record holder is 
PSR J1748— 2446ad, which rotates at 716 Hz lHessels et al.l (120061 ) . This star is tantalizingly 
close to the centrifugal disruption spin frequency for several EOSs (Fig. [1]). Although this 
spin rate can't presently constrain the neutron star radius or the EOS, this situation could 
change in the near future, as more fast-spinning pulsars are discovered. Especially exciting 
would be the discovery of a pulsar near the centrifugal limit whose mass was measured 
through other means; this would lead to limits on the neutron star radius. 



3.2 Mass limits from double- neutron-star systems 

The best EOS constraints to date come from precise measurements of neutron star masses. 
Constraints on masses are obtained by measuring one or more post-Keplerian (PK) (rel- 
ativistic) orbital elements in binary pulsars: advance of periastron, varying redshift and 
time dilation, orbital decay due to gravitational radiation emission, and range and shape of 
Shapiro delay (lTaylorlll992l ). Measurement of two PK parameters uniquely determines the 
masses of the pulsar and companion star. (Measurement of more than two PK parameters 
overconstrains the system and can be used for tests of gravitation; this crucial application 
of binary pulsars is the subject of a separate white paper.) In some sources, most notably 
the double pulsar, additional constraints come from observations of the companion star. 

The most precise mass limits come from double neutron star (DNS) systems, like the 
Hulse- Taylo r binary, PSR B1913-H6 . The pulsar in this system has a mass of (1.4408 ± 
0.0003) M0 (jWeisberg fc Taylorll2003l ). This means that if a EOS predicts a maximum pulsar 
mass below 1.44 Mq, it is automatically excluded. The extremely soft EOS GSl (see Fig. [T]) 
is therefore not a realistic description of matter at supra-nuclear densities. 

The most massive observed neutron star also sets a quantitative upper limit to its central 
density, which then cannot be exceeded in any neutron star. As successively more massive 
neutron stars are observed, the limiting central density gets smaller. Approximately, the 
limit is pmax < 36ps(M0/Mmaa:)^ wherc ps = 2/7 x lO^'' g cm~^ is the nuclear saturation 
density. A large enough observed mass c ould rule out the appeara nce of exotic phases in 
neutron stars, such as deconfined quarks (ILattimer fc Prakashll2005l ). 

Since the discovery of PSR B1913+16 in 1974, at least eight other DNS systems have been 
discovered. The well-measured pulsar masses measured range between 1.24 and 1.44 Mq, a 
surprisingly narrow range, all near the Chandrasekhar limit. 



3.3 Mass limits from other pulsar binaries 

Pulsars in neutron star-neutron star binaries discussed in the previous section have been 
mildly spun up, achieving rotation periods of 10s of milliseconds. By contrast, pulsars in 
neutron star-white dwarf binaries undergo much longer episodes of accretion and are spun 
up to periods of a few milliseconds (hence the name "millisecond pulsar," or MSP). Their 
fast periods lead to superior timing precision. However, they provide difficult targets for PK 
parameter measurements because of very low eccentricities, wide orbits, and low companion 
masses. 

Nevert heless, some measure r nents have been mad e in these systern s. See, for example, 
figure 3 of ILattimer &: PrakashI (120071 ) and figure 3 of iNice et al.l (120081 ). In most cases, the 
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uncertainties on the mass measurements are frustratingly large. In cases with unusually high 
eccentricity or high timing precision, more precise measurements can be made. In the case of 
PSR J0024— 7204H, an MSP in a mildly eccentric (e = 0.07) binary system in globular cluster 
(GC) 47 Tucanae, measurement of periastron advan ce, together with the mass function of 



the system, yields Mp < 1.5 Mq flFreire et all 120031 ). In the case of PSR J1909-3744, a 



bright MSP with a fortuitously "edge -on" geometry, measurement of Shapiro delay yields 
Mp = lA7to°olMQ flHotan et al.ll2006h . These results established that at least some MSPs 



can be recycled with relatively small amounts of mass, but introduced no significantly new 
EOS constraints. 

The situation has now cha nged dramatically. Th e most important event has been the 



discovery of PSR J1903+0327 (jChampion et al.ll2008l ). This extraordinary binary pulsar has 



a period similar to NS-WD systems (P =2.15 ms), but is in a wide {Pf, = 95d), eccentric 
orbit. It was the first MSP discovered in the ongoing ALFA pulsar survey. It is unique 
among fast MSPs in having a fairly massive companion (either a massive white dwarf or 
possibly a main-sequence star, as suggested by the spatial coincidence of a Sun-hke star with 
the astrometric position of the pulsar) and absolutely extraordinary in having an orbital 
eccentricity of 0.44. These properties are completely at odds with all expectations based on 
our present understanding of the stellar evolution in binary systems. 

In addition to being an intriguing system from a binary evolution point of view, the 
unusual orbital characteristics are extremely useful for measuring the pulsar mass. By com- 
bining a precise measurement of the apsidal motion of this system {uj = 86.58 ± 0.25 arc- 
seconds per century) with Shapiro delay measurements, the mass has been determined to 
be 1.671 ± O.OO8M0 (the Aa uncertainty is only O.O32M0). Shapiro delay measurements are 
beginning to over-constrain the system and confirm the unusually high mass of this neutron 
star. The precision of this measurement is unprecedented for any MSP, and the mass is sig- 
nificantly higher than any previous precise neutron star mass measurement. If these results 
stand under further scrutiny, we can exclude several soft EOSs (see red line in Fig. [T]). 

3.4 Millisecond Pulsars in Globular Clusters 

There has been a wealth of MSP discoveries in globular clusters in recent years, thanks to 
deep observations with the Green Bank Telescope and elsewhere; about 2/3 of the MSPs 
now known are located in Because of perturbations from passing stars, or in some 

cases exchange encounters, these systems can be considerably more eccentric than MSPs in 
the Galactic disk, as in the case of PSR J0024— 7204H mentioned above. In fact in some 



cases these MSP binaries are more eccentric than any system in the Galactic disk lFreire et al. 



(I2OO4J . 120071 ). The number of highly eccentric (e > 0.3) MSP binaries increased from 1 known 
before 2004 to 14 at present. Follow-up observations of these pulsars over the coming years 
promise to yield a wealth of PK parameters and, therefore, measurements and constraints 
on neutron star masses. 

Even now, a statistical analysis of these measurements strongly suggests that some 
of these pulsars are massive. By considering the likelihood of different orientations of 
PSRs J1747— 24461 and J in Terzan 5, it has been found that there is a 95% probability 



^see http: //www. naic.edu/ ~pfreire/GCpsr.html| 
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Figure 2: Constraints on the masses of 
PSR J1903+0327 and its companion, based 
on Green Bank and Arecibo timing. The 
hashed area is excluded by knowledge of the 
mass function and by |sini| < 1, where i 
is the orbital inclination. The remaining 
constraints are derived from the rate of ad- 
vance of periastron (u) and the "range" and 
"shape" (r, s) of the Shapiro delay, assuming 
they are of general-relativistic origin. The 
different constraints give a consistent mea- 
surement of the components masses, i.e., GR 
passes this test. The companion is signifi- 
cantly lighter than any NS measured to date, 
it is likely to be a massive WD or a MS star. 
The pulsar is significantly more massive than 
any other NS with a precise mass determina- 
tion. Marginal panels, top: probability dis- 
tribution function for the mass of the pul- 
sar, the median is indicated with a vertical 
line, right: probability distribution function 
for the mass of the companion. From Freire 
et al. (2009). 



that at least one of these objects is more massive than 1.68 M© flRansom et al.l 120051). Fo r 
PSR B1516+02B, in M5, there is a 95% chance that Mp > 1.72 fiFreire et al.ll2008bh . 
Finally, for PSR J1 748-2021B in NGC 6440 there is a 99% probability of Mp > 2.0 



Freire et al.l (l2008al ). These probabilities are all calculated assuming that the apsidal motion 



is purely relativistic. In the case of M5B the non-identification of the companion in HST 
images guarantees that the apsidal motion is relativistic; this is not guaranteed in the case 
of NGC 6440B. A deep HST observation of the environs of NGC 6440B is planned in order 
to investigate the nature of its companion. 

The apsidal motions are now being measured for four more systems in Terzan 5 and two 
others in M28, and it is likely that measurements will be extended to more than a dozen GC 
binaries in the very near future. This, coupled with more mult i- wavelength observations and 
improved statistical methods will allow more detailed knowledge of the underlying distribu- 
tion of MSP masses. Further, we are now reaching the stage at which other PK parameters, 
particularly 7, the time dilation and gravitational redshift parameter, may be measurable in 
some sources, particularly Terzan 5 I. 



3.5 Moment of Inertia of PSR J0737-3039A 

Continued timing of the double pulsar system, J0737— 3039, will greatly improve the pre- 
cision of its PK parameters, in particular the measurement of the orbital decay due to the 
emission of gravitational waves. Eventually, the precision of this measurement is limited by 
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uncertainties in the kinetic contribution to the orbital decay. Present (Deller, Bailes and 
Tingay, larXiv:0 902.0996) and future interferometric studies will allow a very precise esti- 
mate of the contribution of kinetic effects to its orbital decay. By doing this we might be 
able to measure the expected spin-orbit contribution to the uj (jPamour fc Schaferl 119881) of 
PSR J0737— 3039A, allowing an estimate of this pulsar's moment of inertia (IKramer &: Wex 



20091 ). This, together with its well-determined mass, will introduce unique ari d fundamental 



constraints on the equation of state for dense matter (ILattimer fc Schutzll2005l ). These would 



complement the constraints derived from the measurement of MSP masses. 



4 Key Requirements 

High Sensitivity Facilities Pulsars are low luminosity sources, and most of the objects 
that are providing the new mass limits, the GC pulsars, are exceptionally faint. Detection 
of PK parameters typically requires consistent observation of a pulsar over a period of 
years. Obtaining any new results on these sources absolutely requires the world's most 
sensitive radio telescopes, particularly the GET and Arecibo Observatory. These are the 
only instruments capable of even detecting the vast majority of the GC pulsars — and the 
only ones capable of detecting PSR J1903-f0327 with the S/N required to measured the 
three PK parameters we now measure for that system. 

Even Higher Sensitivity Facilities The additional sensitivity of an instrument like the 
Square Kilometer Array (SKA) would revolutionize this field, by vastly increasing the 
sample of MSPs with measurable masses and greatly increasing the number of measurable 
PK parameters in all sorts of systems. It will also be able to determine the distance to 
PSR J0737— 3039 with enough precision to allow a precise measurement of the moment of 
inertia of PSR J0737-3039A. 

Improved Radio Instrumentation Precision pulsar timing requires observing the pulsar 
signal with high time resolution over large radio bandwidths (~1 GHz). The development 
of the GET and upgrades to Arecibo, have made such wide bandwidths available for 
the first time, and back-end instrumentation is continually evolving toward taking full 
advantage of these bandwidths. This increase in bandwidth and time resolution is essential 
to fully resolve the radio pulses' fine structure and increase the sensitivity to these faint 
sources. 

Continuous Long-term Observations For all GC pulsars, only continued, long-term 
timing will allow measurement of the relativistic 7, which is vital for precise measure- 
ments of pulsar masses. In some important cases, like NGC 6440E, this will require a 
couple of decades. Improvements in sensitivity and timing precision will be important 
to reduce these time horizons, but the long-term nature of these projects is essentially 
inevitable. 

New MSPs Progress in using pulsar mass measurements to constrain the EOS of neutron 
star matter has only been possible because of the discovery of new objects, like PSR 
J1903+0327 or the GC pulsars. This is driven by continual improvement in telescope 
sensitivity, available bandwidth, computing capabilities, and pulsar search algorithms. 
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Continued searches, like the ALFA pulsar survey, the 350-MHz drift-scan survey with the 
GBT, the 327-MHz drift-scan survey with Arccibo and others might discover new binary 
MSPs with orbital characteristics that make their masses easy to measure. Some of them 
might have masses that happen to constrain the EOS. Whatever their masses, the more we 
increase the sample of MSPs with measurable masses, the better will be our understanding 
of their mass distribution. 

The improvements in time resolution that allow better timing precision also allow much 
higher search sensitivity to fast-spinning pulsars. The present pulsar surveys can detect 
sub-millisecond pulsars in Galactic volumes that arc hundreds of times larger than any 
previous pulsar surveys. As discussed above, the discovery of such systems also has the 
potential to fundamentally constrain the EOS. 
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